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Isovector EMC effect explains the NuTeV anomaly 
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A neutron or proton excess in nuclei leads to an isovector-vector mean-field which, through its 
coupling to the quarks in a bound nucleon, implies a shift in the quark distributions with respect 
to the Bjorken scaling variable. We show that this result leads to an additional correction to the 
NuTeV measurement of sin'^ Qw- The sign of this correction is largely model independent and acts 
to reduce their result. Explicit calculation within a covariant and confining Nambu-Jona-Lasinio 
model predicts that this vector field correction accounts for approximately two-thirds of the NuTeV 
anomaly. We are therefore led to offer a new interpretation of the NuTeV measurement, namely, 
that it is further evidence for the medium modification of the bound nucleon wavefunction. 

PACS numbers: 24.85.-|-p, 13.60.Hb, ll.80.Jy, 21.65.Cd 



Within relativistic, quark-level models of nuclear struc- 
ture, the mean scalar and vector fields in the medium 
generate fundamental changes in the internal structure 
of bound hadrons. These modifications lead to a good 
description of the EMC effect in finite nuclei and predict 
a more dramatic modification of the bound nucleon spin 
structure function [1-3]. We show that in nuclei with 
N Z this approach leads to interesting and hitherto 
unexplored effects connected with the isovector-vector 
mean- field, which is usually represented by the , and is 
in part responsible for the symmetry energy. In a nucleus 
like ^^Fe or ^ospb where N > Z, the p° field will cause 
the M-quark to feel a small additional vector attraction 
and the d-quark to feel additional repulsion. 

In this Letter we explore the way in which this addi- 
tional vector field modifies the traditional EMC effect. 
However, there is an even more important issue which is 
our main focus. Even though the mean-field is com- 
pletely consistent with charge symmetry, the familiar as- 
sumption that Up{x) = dn(x) and dp(x) = Un{x) will 
clearly fail for a nucleon bound in a nucleus with N Z. 
Therefore correcting for the p° field is absolutely critical 
in a situation where symmetry arguments are essential, 
such as the use oi N ^ Z nuclear data from v and V 
deep inelastic scattering (DIS) to extract sin^ Qw via the 
Paschos-Wolfenstein relation [4]. Indeed, we show that 
the deviation from the naive application of charge sym- 
metry to the v and V data on ^^Fe naturally explains the 
famous NuTeV anomaly. 

The Paschos-Wolfenstein (PW) ratio is defined by [5] 
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where A represents the target, NC indicates weak neu- 
tral current and CC weak charged current interaction. 
Expressing the cross-sections in terms of quark distribu- 
tions and ignoring heavy flavour contributions, the PW 



ratio becomes 
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where xa is the Bjorken scaling variable of the nucleus 
multiplied by A, (. . .) implies integration over xa, and 
= qA — qA are the non-singlet quark distributions of 
the target. Therefore, the quantities in the angle brackets 
are simply the momentum fractions of the target carried 
by the valence quarks. 

Ignoring quark mass differences and possible elec- 
troweak corrections the u- and d-quark distributions of an 
isoscalar target will be identical, and in this limit Eq. (2) 
becomes 
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If corrections to Eq. (3) are small the PW ratio provides 
a unique way to measure the Weinberg angle. 

Motivated by Eq. (3) the NuTeV collaboration ex- 
tracted a value of sin^ Qw from neutrino and anti- 
neutrino DIS on an iron target [6], finding 

sin^ Qw = 0.2277 ± 0.0013(stat.) ± 0.0009(syst.). (4) 

The three-sigma discrepancy between this result and the 
world average [7], namely sin^ Qw = 0.2227 ± 0.0004, is 
the NuTeV anomaly. Some authors have speculated that 
the NuTeV anomaly supports the existence of physics 
beyond the Standard Model [8]. However, existing high 
precision data for other electroweak observables places 
tight constraints on new physics explanations. Standard 
Model corrections to the NuTeV result have largely been 
focused on nucleon charge symmetry violating effects [9] 
and a non-perturbative strange quark sea [8]. Charge 
symmetry violation, arising from the u- and d-quark mass 
differences, is probably the best understood and con- 
strained correction and can explain approximately one- 
third of the NuTeV anomaly [10]. Standard nuclear cor- 
rections like Fermi motion and binding are found to be 
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small [11] . However effects from the medium modification 
of the bound nucleon, which are now widely accepted as 
an essential ingredient in explaining the EMC effect [12], 
have hitherto not been explored in relation to the NuTeV 
anomaly. 

In our approach, presented in Refs. [2, 3, 13], the scalar 
and vector mcan-ficlds inside a nucleus couple to the 
quarks in the bound nucleons and self-consistently mod- 
ify their internal structure. The scalar field renormalizes 
the constituent quark mass, resulting in effective hadron 
masses in-medium. The influence of the vector fields on 
the quark distributions arises from the non-local nature 
of the quark bilinear in their definition [13]. This leads 
to a largely model independent result for the modifica- 
tion of the in-medium parton distributions of a bound 
nucleon by the vector mean-fields [13-15], namely 
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The subscript indicates the absence of vector fields and 
p+ is the nucleon lightcone plus component of momen- 
tum. The quantities and are the lightcone plus 
component of the net vector field felt by the nucleon and 
a quark of flavour q, respectively. Importantly Eq. (5) is 
consistent with baryon number and momentum conser- 
vation, and implies that the mean vector field carries no 
momentum. 

Before embarking on explicit calculations, we first ex- 
plore the model independent consequences of Eq. (5) for 
the PW ratio and the subsequent NuTeV measurement 
of sin^ Ow The NuTeV experiment was performed on 
a predominately ^^Fe target, and therefore isoscalarity 
corrections need to be applied to the PW ratio before 
extracting sin^Ov^. Isoscalarity corrections to Eq. (3) 
for small isospin asymmetry have the general form 
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where the dependence of this correction resides com- 
pletely with sin^ Qw NuTeV perform what we term 
naive isoscalarity corrections, where the neutron excess 
correction is determined by assuming the target is com- 
posed of free nucleons [16]. However there are also 
isoscalarity corrections from medium effects, in particular 
from the medium modification of the structure functions 
of every nucleon in the nucleus, arising from the isovec- 
tor p° field. For nuclei with N > Z the p° field develops 
a non-zero expectation value that results in Vu < Vd, 
that is, the w-quarks feel less vector repulsion than the 
d-quarks. A direct consequence of this and the transfor- 
mation given in Eq. (5) is that there must be a small shift 
in quark momentum from the u- to the d-quarks. There- 
fore the momentum fraction {xa — a;^ rf^) in Eq. (6) 
will be negative, even after naive isoscalarity corrections 
are applied. Correcting for the p'^ field will therefore have 



the model independent effect of reducing the NuTeV re- 
sult for sin^ Qw- As we shall see, this correction largely 
explains the NuTeV anomaly. 

To determine the nuclear quark distributions we use 
the Nambu-Jona-Lasinio (NJL) model [17], which is 
viewed as a low energy chiral effective theory of QCD 
and is characterized by a 4-fermion contact interaction 
between the quarks. The NJL model has a long his- 
tory of success in describing mesons as qq bound states 
[18] and more recently as a self-consistent model for free 
and in- medium baryons [2, 3, 13, 19]. The original 4- 
fermion interaction term in the NJL Lagrangian can be 
decomposed into various qq and qq interaction channels 
via Fierz transformations [20] . The relevant terms of the 
NJL Lagrangian to this discussion are 
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where Z?-^ = yfA^^ {A = 2,5,7) are the the colour 3 

matrices [2] , (7 = ^7270 and m is the current quark mass. 

The scalar qq interaction term generates the scalar 
field, which dynamically generates a constituent quark 
mass via the gap equation. The vector qq interaction 
terms arc used to generate the isoscalar- vector, cjq, and 
iso vector- vector, poj mean-fields in- medium. The qq in- 
teraction terms give the diquark f-matrices whose poles 
correspond to the scalar and axial- vector diquark masses. 
The nucleon vertex function and mass are obtained by 
solving the homogeneous Faddeev equation for a quark 
and a diquark, where the static approximation is used to 
truncate the quark exchange kernel [19]. To regularize 
the NJL model we choose the proper-time scheme, which 
enables the removal of unphysical thresholds for nucleon 
decay into quarks, and hence simulates an important as- 
pect of confinement [21, 22]. 

To self-consistently determine the strength of the mean 
scalar and vector fields, an equation of state for nuclear 
matter is derived from the NJL Lagrangian, Eq. (7), us- 
ing hadronization techniques [22]. In a mean-field ap- 
proximation the result for the energy density is [22] 
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where the vacuum energy Ey has the familiar Mexican 
hat shape and the energies of the protons and neutrons 
moving through the mean scalar and vector fields are 
labelled by £p and f„, respectively. The corresponding 
proton and neutron Fermi energies are 



epa = EFa+ Va = ^ + p%^ + 3UJ0 ± po, (9) 
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FIG. 1: Isospin dependence of the EMC effect for proton- 
neutron ratios greater than one. The data is from Ref. [24] 
and corresponds to N = Z nuclear matter. 



where a = p or n, the plus sign refers to the proton, 
is the in- medium nucleon mass and ppa the nucleon 
Fermi momentum. Minimizing the effective potential 
with respect to each vector field gives the following useful 
relations: uq = 6Guj [pp + Pn) and po = 2Gp {pp ~ pn), 
where pp is the proton and pn the neutron density. The 
vector field experienced by each quark flavour is given by 
K = t^o + Po and Vd = ujq - Po- 

As explained in Ref. [2] , the parameters of the model 
are determined by standard hadronic properties, and the 
empirical saturation energy and density of symmetric nu- 
clear matter. The new feature of this work is the 
field, where Gp is determined by the empirical symmetry 
energy of nuclear matter, namely 04 = 32 MeV, giving 
Gp = 14.2GeV-2. 

Details of our results for the free and N ~ Z in- medium 
parton distributions are given in Refs. [2, 3, 19]. For in- 
medium isospin dependent parton distributions our pro- 
duce is as follows: Effects from the scalar mean-field are 
included by replacing the free masses with the effective 
masses in the expressions for the free parton distributions 
discussed in Ref. [19]. For asymmetric nuclear matter 
[N ^ Z) the proton and neutron Fermi momentum will 
differ and therefore so will their Fermi smearing func- 
tions. To include the nucleon Fermi motion, the quark 
distributions modified by the scalar field are convoluted 
with the appropriate Fermi smearing function, namely 
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where J\fp = Z. 
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Vector field effects can be included in Eq. (10) by the sub- 
stitutions Epa ^ £Fa and Mn -^Mn = ^SFp + ^^Fn- 
Our final result for the infinite asymmetric nuclear mat- 
ter quark distributions, which includes vector field effects 
on both the quark distributions in the bound nucleon and 



FIG. 2: Isospin dependence of the EMC effect for proton- 
neutron ratios less than one. 

on the nucleon smearing functions, is given by 
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The subscript ^0 indicates a distribution which includes 
effects from Fermi motion and the scalar mean-field. The 
distributions calculated in this way are then evolved [23] 
from the model scale, Qg = 0.16 GeV^, to an appropriate 
for comparison with experimental data. 
The EMC effect is defined by the ratio 
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where qa are the quark distributions of the target and qf 
are the distributions of the target if it was composed of 
free nucleons. Results for the isospin dependence of the 
EMC effect are given in Figs. 1 and 2. 

Fig. 1 illustrates the EMC effect for proton rich matter, 
where we find a decreasing effect as Z/N increases. An 
intuitive understanding of this result may be obtained by 
realizing that it is a consequence of binding effects at the 
quark level. For Z/N > 1 the po field is positive, which 
means Vu > Vd and hence the w-quarks are less bound 
than the d-quarks. Therefore the u-quark distribution 
becomes less modified while medium modification of the 
d-quark distribution is enhanced. Since the EMC effect 
is dominated by the it-quarks it decreases. The isospin 
dependence of the EMC effect for nuclear matter with 
Z/N < 1 is given in Fig. 2. Here the medium modifica- 
tion of the u-quark distribution is enhanced, while the d- 
quark distribution is modified less by the medium. Since 
the EMC ratio is initially dominated by the w-quarks the 
EMC effect first increases as Z/N decreases from one. 
However, eventually the d-quark distribution dominates 
the ratio and at this stage the EMC effect begins to de- 
crease in the valence quark region. We find a maximal 
EMC effect for Z/N ~ 0.6, which is slightly less than the 
proton-neutron ratio in lead. This isospin dependence is 
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clearly an important factor in understanding the A de- 
pendence of the EMC effect, even after standard neutron 
excess corrections are applied. 

Now wc turn to the consequences of the isospin depen- 
dence of the EMC effect for the NuTeV measurement of 
sin^O^y. The NuTeV experiment was performed on an 
iron target, which, because of impurities had a neutron 
excess of 5.74% [6]. Choosing our Z/N ratio to give the 
same neutron excess, we use our medium modified quark 
distributions and Eq. (6) to determine the full isoscalar- 
ity correction to the isoscalar PW ratio, given by Eq. (3). 
Using the Standard Model value for the Weinberg angle 
we obtain ARpw = —0.0139. If we break this result 
into the three separate isoscalarity corrections, by using 
Eq. (6) and the various stages of the medium modified 
quark distributions, we find 

ARpw = AR'^^'^ + Ai?p™' -I- ARp^ 

= -(0.0107 + 0.0004-1- 0.0028). (13) 

Higher order corrections to Eq. (6) do not change this 
result. The NuTeV analysis includes the naive isoscalar- 
ity correction [16] but is missing the medium correction 
of —0.0032 [25]. This new correction accounts for two- 
thirds of the NuTeV anomaly. If we also include the 
well constrained charge symmetry violation (CSV) cor- 
rection, AR$I^ = -0.0017 [10], which originates from 
the quark mass differences, we have a total correction of 
AiiSS^i"™ + AR9,-^ = -0.0049. The combined correc- 



tion accounts for the NuTeV anomaly [27]. 

Since CSV and medium modification corrections ex- 
plain the discrepancy between the NuTeV result and the 
Standard Model, we propose that this NuTeV measure- 
ment provides strong evidence that the nucleon is mod- 
ified by the nuclear medium, and should not be inter- 
preted as an indication of physics beyond the Standard 
Model. In our opinion this conclusion is equally profound 
since it may have fundamental consequences for our un- 
derstanding of traditional nuclear physics. We stress that 
the physics presented in this paper, in particular the ef- 
fects of the p'^ mean-field, are consistent with existing 
data [28], but can strongly influence other observables. 
For example, the field gives rise to a strong flavour 
dependence of the EMC effect, and experimental propos- 
als have been submitted at Jefferson Lab to look for such 
eff'ects. 
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